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Abstract 

Using an improved estimator in the loop-cluster algorithm, we inves- 
tigate the constraint effective potential of the magnetization in the spin 
~ quantum XY model. The numerical results are in excellent agreement 
with the predictions of the corresponding low-energy effective field the- 
ory. After its low-energy parameters have been determined with better 
than permille precision, the effective theory makes accurate predictions 
for the constraint effective potential which are in excellent agreement 
with the Monte Carlo data. This shows that the effective theory indeed 
describes the physics in the low-energy regime quantitatively correctly. 
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1 Introduction 



When a nonperturbative system of many strongly coupled degrees of freedom under- 
goes the spontaneous breakdown of a continuous global symmetry, massless Gold- 
stone bosons arise as the relevant low-energy degrees of freedom. Despite the fact 
that the underlying microscopic system can usually only be studied numerically, the 
low-energy dynamics of the Goldstone bosons can be addressed analytically using 
a systematic low-energy effective field theory. The a priori unknown low-energy 
parameters of the effective theory can then be determined by matching the results 
of numerical simulations of the underlying microscopic system to analytic results of 
the effective field theory. Thanks to this interplay between numerical simulations 
and analytic effective field theory calculations, important insights have been gained 
into the Higgs sector of the Standard model [H, 2[ and the dynamics of QCD as 
well as of magnetic systems, including the undoped precursors of high-temperature 
superconductors [!, 0]. The latter are described by a low-energy effective theory 
for magnons — the Goldstone bosons of the spontaneously broken SU (2) spin 
symmetry. Using this theory, the shape of the constraint effective potential of the 



order parameter has been worked out in detail by Gockeler and Leutwyler [10|, 111 
Their predictions for the shape of the constraint effective potential had already 
been tested against Monte Carlo simulations of classical 3-d 0(3) and 4-d 0(4) lat- 



tice models [12l . Recently, we have performed a high-accuracy investigation of the 



antiferromagnetic spin | quantum Heisenberg model on a square lattice [13j. In 
particular, we have employed a new improved estimator for the distribution of the 
staggered magnetization using the loop-cluster algorithm 0, 0, Q3]. The very accu- 
rate Monte Carlo data were then compared with the analytic results of the magnon 
effective theory, resulting in a determination of the low-energy parameters with per- 
mille accuracy. In this paper, we extend this investigation to the (2 + l)-d quantum 
XY model, which models quantum magnetism, as well as superfluidity of hard-core 
bosons and has a spontaneously broken £7(1) symmetry. The analytic results of 
the effective field theory readily apply to the Goldstone bosons of superfluidity. In 
the present paper, however, we will use "magnetic" language throughout, having in 
mind quantum magnetism. After its low-energy parameters have been determined 
with better than permille precision, the effective theory makes unambiguous pre- 
dictions for the constraint effective potential, which are in excellent agreement with 
the Monte Carlo data. 

The rest of this paper is organized as follows. In section 2 we summarize the 
predictions of the low-energy effective field theory. In section 3 we present the 
results of our numerical simulations obtained with the loop-cluster algorithm, which 
are compared with the effective field theory predictions in section 4. Finally, section 
5 contains our conclusions. 
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2 Effective Field Theory Predictions 



In section 4 we will compare our very accurate Monte Carlo data with the effective 



field theory predictions of Gockeler and Leutwyler [10 . Illj which are derived from a 
scalar 3-d 0(A r )-symmetric effective field theory. Here, we summarize those results 
that are relevant for our study. While these results were derived in the framework 
of a 3-d relativistic quantum field theory, they are readily applicable to the (2 + 1)- 
d quantum XY model, because at low energies its Goldstone boson has a linear 
"relativistic" dispersion relation. The Hamiltonian of the quantum XY model is 
defined as 

H = -JJ2( S 1 S 1 + S l S D -M-B, (2.1) 

i x v) 

where x and y are nearest-neighbor sites on a square lattice with spacing a, and 
J > is a constant. Working in natural units in which h—1, the spin | operators 
S x obey the standard commutation relations 

[S*,S b y ]=i5 xy e abc S c x . (2.2) 

The quantity B = (Bi, Bq) is a uniform magnetic field in the XY plane that couples 
to the magnetization order parameter 



M 



(£^'£<)- (2-3) 



In the infinite volume limit and at zero temperature, the vacuum expectation value 
of M is non-zero, signaling the spontaneous breakdown of the U(l) spin symmetry, 
which is generated by the 3-component of the total spin 



S 3 = J2 # ( 2 - 4 ) 



It should be noted that, in contrast to the Heisenberg model which has an SU(2) 
spin symmetry, even at B — 0, here only the generator S 3 commutes with the 
Hamiltonian. It should also be pointed out that, on a bipartite lattice, in the 
absence of a magnetic field, both the antiferromagnetic and the ferromagnetic XY 
model describe the same physics, since they are related to each other by a unitary 
transformation. 

At low energies, the relevant degrees of freedom are the Goldstone bosons result- 
ing from the spontaneously broken global symmetry. In the present case of the XY 
model, the U(l) spin symmetry breaks completely, and we thus have one magnon (or 
superfluid Goldstone boson, depending on the physical interpretation of the quan- 
tum XY model). The low-energy effective field theory is formulated in terms of the 
magnetization order parameter field 

e(x) = ( ei (x), e 2 (x)) G S\ e(x) 2 = 1, (2.5) 
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where x = (xi,X2,t) is a point in Euclidean space-time. Up to higher-order correc- 
tions in a derivative expansion, the effective action takes the form 



S[e\ = / d 2 x dt 



| (die ■ d t e + ^d t e- d t e ) - Mr ■ D 



(2.6) 



where p is the spin stiffness, c is the spin-wave velocity, and M is the magnetization 
density. The partition function is then given by 



J De exp(-S[e[). (2.7) 



The dispersion relation of the magnon takes a "relativistic" form, where the velocity 
of light is replaced by the spin-wave velocity c. By introducing £3 = ct, the effective 
action can thus be written in the manifestly Euclidean space-time rotation-invariant 
form 

5[e] = J d 3 x - c (|d M e • <9 M e - Me ■ B) . (2.8) 

Since the underlying XY model does not exhibit Euclidean rotational invariance, 
the symmetry is accidental and exists only in the leading terms of the derivative 
expansion. Indeed, higher-order four-derivative terms in the effective action break 
Euclidean space-time rotation invariance down to the discrete rotation subgroup of 
the square lattice. 

The system, described by the effective action ( 12. 6p . is considered in a periodic 
cubic space-time volume L x L x j3 with the inverse temperature fixed at (3 = L/c^ 
The space-time average of the magnetization is given by 

$ = \h& I d2x dt e>) = \hj d ' x e>) - (2,9) 



In contrast to llOl . 1 1 11 1 , we have included a factor ~ in the definition of $ because 



the quantum spins of the underlying XY model have S — |, while the effective field 
e(x) is normalized to 1. The probability distribution of the mean magnetization $ 
is obtained as a 5-function constrained path integral for the partition function 

= \ j De exp(-S[e1) 5 ($ -~ J d 3 x %x)) . (2.10) 

As a consequence of the U(l) symmetry, it only depends on the magnitude $ = |$|. 
This distribution is normalized according to 

POO 

2tt / rf$ $ p($) = 1. (2.11) 







One of the fundamental quantities in the present study is the constraint effective 
potential «($). It represents the free energy density of the model which is obtained 



1 We set the Boltzmann constant ks = 1. 
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by computing the path integral over configurations constrained to a given fixed mean 
magnetization value $, and is determined by 



= A/"exp(-L 3 u($)). 



(2.12) 



The analytical expression for the normalization factor A/ - , derived in lOj reads 



pL 



M2 47r 2 ce /3 °/ 2 



1 + 



L 2 



(2.13) 



where we have defined the magnetization per spin Ai = Aia 2 . The quantity /3q is a 
shape- dependent coefficient characterizing the geometry of the space-time box. For 
the exactly cubical space-time volume considered here it is given by 0o = 1.45385. 
In the infinite-volume and zero-temperature limit the constraint effective potential 
approaches the infinite volume effective potential which is known to be a convex 
function of $ 15, [l6j]. In a finite volume, on the other hand, w($) is not necessarily 



convex. We may define an extensive variant of the intensive quantity «($) as 

U(<$>) = L 3 w($). (2.14) 



Within the effective field theory framework, the finite-size corrections to the con- 
straint effective potential were systematically worked out by Gockeler and Leutwyler 



10l . Illfl . Near its minimum the 1/L expansion of t/($) takes the form 



U{*) = U (ij) + -jL UM + O (±\ . (2.15) 

The quantities Uo(ip) and Ui{ip) depend on L only through the rescaled variable 

pL $ - M , 

i) = - 2.16 

c M 

The leading order contribution to the constraint effective potential is given by the 
inverse Laplace transform 



exp(— £/oC0)) = / dx exp(—ixvp + T(ix)) (2-17) 

of the function 



n=0 



Again, the quantities (3 n are shape- dependent coefficients of the finite space-time box 
which are described in detail in appendix B of [l7[. Remarkably, the function T(ix) 
is entirely kinematical and thus, unlike Ui(ip), the quantity U (ip) is universal, i.e. 
completely independent of the low-energy parameters 0. Hence U (?P) is the same 
for all 3-d systems with a spontaneously broken U(l) symmetry, including the 3-d 
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classical and the (2 + l)-d quantum XY model. The 1/L correction to the leading 
contribution U (ip) is given by 

POO 

Uity) = if) + exp([/ (») / dxRe{exp(-ixip + T(ix))£l(ix)}, (2.19) 

Jo 

with 

m = -~ (W - mo - 1^3) - ^ 2 , 

00 „ 

w (0 = E^3iv^" 1 - (2 ' 20) 

n=l V 

Furthermore, fco appearing above is a low-energy constant which is given by 

* = + + (2 ' 21) 

where /ii and /12 are the low-energy constants associated with the higher-order terms 
in the effective action 



AS[e\ = - J d 2 x dt ht{e- B) 2 + h 2 B 2 . (2.22 



A non-zero magnetic field B in the XY plane turns the magnons into pseudo-Nambu- 
Goldstone bosons with a non-zero mass m determined at leading order by 



MB 



m 2 = — -, B = \B\. (2.23) 



The constant ko also appears in the 5-dependence of the field expectation value 

8\oLJ ^ nAnol 



8\pL/ ' nAn 2 \ 

r ni,n 2 =0 1 - 

^ n=0 



v 2n-2 



1 C ^ /3„+l , , 2 n 1/ C 2 1 

> — —(mcL) — 

2pL^ n! 1 J 8>L ; ( mc L) 

1 c 1 1 fmc 2 \ 2 fmc 
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-6^) 2 W^)W) • (2,4) 



2pL(mcLy 64tt 2 V p 

It should be noted that eq. fl2.24p was derived in the p- regime of chiral perturbation 
theory in which mcL 3> 1 while mc 2 , c/L <C 27rp. In particular, in eq.f )2.24p one 
cannot make B (and thus m) arbitrarily small, because one would otherwise enter 
the e- regime in which mcL w 1. 
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The low-energy constant k can be determined either from Ui(ip) or by 
fitting the relevant Monte Carlo data to the corresponding theoretical predictions 
(eq. fl2.19p and eq. fl2.24p ). As we will demonstrate later, the numerical values for k 
obtained from Ui(ip) and \ (§)(B)\ are consistent. 

Gockeler and Leutwyler have also worked out analytic predictions for the first 
and second moment of the probability distribution p($) up to two loops. They 
obtained 

m = M( l + ±?i + jLS\ + o(l. 

(C-W)")^!^^). P.25) 

where the additional shape- dependent coefficients for the cubic box considered here 
are given by (3 l = 0.225785 and (3 2 = 0.0 1 06 08 [l7 . 

Other physical quantities of central interest are the susceptibilities. First, one 
identifies the order parameter susceptibility 

Xx = j2[ dt |Tr[M 1 (0)M 1 (t)exp(-/3i/)]. (2.26) 

Here M 1 = is the first component of the magnetization. A second suscepti- 

bility refers to the U{1) conserved quantity M 3 and is defined as 



X3 



l?L dt ^ 1& [ M3 ( ) Af3 (*) ex P(-^)]> ( 2 - 27 ) 

with M 3 = J2 X Both xi an d X3 can be measured very efficiently with the 
loop-cluster algorithm using improved estimators jsj. 

Another reference that provides analytic effective field theory results, which can 
be compared with our Monte Carlo data, is a paper by Hasenfratz and Niedermayer 
(9)]. Using magnon chiral perturbation theory up to two-loop order, they obtained 
the finite-size and finite-temperature effects of Xi i n the e-regime 

Xi = { 1 + -T7/M0 + U^l) 2 NO 2 + A(0] + O (±) ) . (2.28) 



pLl v ' 2 \pLl 

The quantity I = ((3c/ L) 1 ^ 3 determines the shape of an approximately cubic space- 
time box of volume L x L x (3, with (3c ~ L. The functions /?»(/) are known 
shape- dependent coefficients [H, 17]. For an exactly cubical space-time volume (i.e. 
for I — 1) the result of eq. fl2.28p agrees with eq. (12.251) since 

(($ - <$» 2 > + <$> 2 = <$ 2 > = §^ = (2.29) 
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The factor 2 arises due to the two components of the magnetization vector. Re- 
markably, up to two-loop order the analogous expression for % 3 takes the simple 
form 

which does not display any corrections of lower orders, neither of 0(1/ L) nor 
0(l/L 2 ). 

The above expressions have been used to determine the low-energy parameters 



by a fit of Xi an d Xz t° Monte Carlo data [18 

M = 0.43561(l)/a 2 , p = 0.26974(5) J, c = 1.1347(2) Ja. (2.31) 

In this very accurate study, the cubical geometry has been reached by tuning j3 
until temporal and spatial winding numbers agreed. The spin-wave velocity has 
then been determined as c = L//3 with fraction of a permille precision. The fitted 



magnetization density is consistent with the result M. = 0.437(2) / a 2 obtained in 19 
For the 2-d spin | Heisenberg model, using the same method, the corresponding low- 
energy parameters have recently also been determined with fraction of a permille 
accuracy in 



3 Probability Distribution of the Magnetization 

The quantum XY model can be simulated very efficiently with the loop-cluster al- 
gorithm 0,0) 14]. In 21] the Kosterlitz-Thouless phase transition has been studied 
for the first time using the loop-cluster algorithm. The transition temperature Tkt 
has been determined very precisely from the winding numbers, which was not pos- 



sible before. By using the same improved estimator as introduced in |13J we extract 
the probability distribution of the magnetization. Every cluster contributes addi- 
tively to the first component of the magnetization. The cluster size |C| (i.e. the 
number of lattice points in a given cluster) determines the first component of the 
magnetization of the cluster C, which is proportional to ±|C|. Under cluster flip 
the magnetization of a cluster changes sign. Starting from a given spin configu- 
ration, the distribution of the magnetization is recorded as a histogram which is 
built iteratively using one cluster after another. The initial histogram pi(m) (with 
m G {— M, —M + l, 0, M — 1, M}, where M is the number of space-time lattice 
points) is constructed from the first cluster as 

Pi( m ) = g [<Wil + S ™-\Ci\] • 

The two entries of the initial histogram correspond to the two possible orientations 
of the first cluster, each arising with probability |. In the i-th iteration step (with 
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Figure 1: Examples of histograms pN^m) obtained for two individual spin configura- 
tions on a 16 2 lattice using the improved estimator. 

i G {2, 3, N}), where N is the number of clusters of a given configuration, a new 
histogram Pi(m) is built from the previous one as 



Pi( m ) = o [Pi-i( m + \ c i\) +Pi-i(m - \Ci\ 



(3.2) 



The final histogram after N steps is given by pj^{m). In figure 1 we show examples of 
histograms ]?/v(m) obtained for two individual spin configurations. In the left panel 
the example contains one cluster that is bigger than all the other clusters together. 
Therefore the region around m = is not sampled. Additionally, there are two 
relatively large clusters that give rise to the multiple peaks in the distribution. On 
the other hand, in the example shown in the right panel, there are two clusters of 
similar size, such that the region around m = is also sampled. The average of 
the histograms pjv(m) for all configurations in the Markov chain generated by the 
cluster algorithm yields the final probability distribution of the first component of 
the magnetization 

p(m) = (pjv(m)). (3.3) 



By construction, it is properly normalized as 



M 



(3.4) 



m=-M 



The numerical effort to build the improved estimator is proportional to the number 
of lattice points M and, in addition, proportional to the number of clusters N. Since 
the number of clusters is proportional to the volume, the evaluation of the improved 
estimator requires a computational effort proportional to M 2 , and thus becomes 
rather time-consuming for large volumes. Still, the improved estimator increases the 
statistics by a factor of 2 N which is exponential in the volume. Therefore investing 
a polynomial effort M 2 should still be justified. Compared to 13j where the method 
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had been implemented in a straightforward manner, we have been able to achieve 
substantial improvements. This allows us to investigate volumes as large as 64 2 
while our original study of the Heisenberg model was limited to 24 2 lattices. 

Remarkably, the computational effort to build the histograms can be reduced 
by a factor of order up to 1000 (for a volume 64 2 ) by several optimizations. First 
of all, it is obvious that one should evaluate eq.(3.2) only for m G {— Mj, —Mi + 
1, 0, Mi -1, MJ with 

M l = J2\Cj\ (3-5) 

for the i-th iteration step. Therefore most of the zero values at the border of the 
partial histograms are not evaluated explicitly. By sorting the clusters by their sizes 
and building the histogram by starting with the smallest clusters and finishing with 
the largest, M, and therefore the computational effort per iteration grows as slowly as 
possible. Due to the time-discretization [3] (which has four Trotter steps), all cluster 
sizes \Ci\ are multiples of 2. Hence the cluster sizes can be divided by 2 when building 
the histograms. This halves memory usage and it also halves the loop sizes for 
evaluating equation (3.2). A further optimization is to use two arrays for the partial 
histograms. One is the source and the other the destination for the new histogram. 
The roles alternate per iteration. This prevents some memory allocations, copy 
operations and initializations with zeros. The histogram is symmetric. So only the 
part for m G {0, M — 1, M} is held in memory and equations (3.1) and (3.2) are 
modified as 

Pi( m ) = 2<Wil> ( 3 - 6 ) 

and 

Pi (m) = \ \pi-i{\m + \d\\) +Pi-i{\m - \&\\)] ■ (3.7) 

Generating a cluster size histogram for one configuration one observes a relatively 
large number of clusters with the same small cluster size. Hence, an improvement 
is to treat clusters of equal sizes in one iteration using the fact, that a histogram of 
an even number of clusters of only one equal size equals to 

n \c\ 

PevenM = Yl \ 2i£i. + fcj *M.2*|C| , (3.8) 

where n\c\ is the number of clusters of size \C\. For odd n\c\ the corresponding 
histogram equals to 



Podd{m 

k=0 



* i c 

(3.9) 

i—n V 9 r KJ 



A further improvement uses the fact that partial histograms built by eq.(3.1) always 
contain zero values either for odd or for even m. This holds for arbitrary cluster 
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sizes \Ci\. By tracking these two cases and using only the non-zero values one can 
further optimize the computational effort. All these optimizations do not influence 
the resulting histogram. For large volumes the computational effort is still growing 
with the square of the volume. 

The final optimization is more delicate. We have found that it is possible to 
divide the cluster size by some volume-dependent factor larger than 2. The round 
off error is treated by an error propagation technique, which is currently not yet 
fully optimized. The division of the cluster size alters the resulting histogram. 
The dividing factor is chosen empirically, such that the resulting systematic error 
is smaller than the statistical error of the Monte Carlo data. For our simulations 
the dividing factor was proportional to L. The computational effort grew with a 
power of about 1.3 of the volume. This method could also be used for simulations 
in continuous time which result in non-integer cluster sizes. 

The mean value of the first component of the magnetization $x corresponding 
to a given value of m is 

777 

$ l = 771T7- (3-10) 

The factor 2 arises because we are dealing with quantum spins |. Now one can iden- 
tify the probability distribution of the first component of the mean magnetization 

as 

p($i) d$! = p{m). (3.11) 

It turns out that the non-zero entries of the histogram p(m) correspond to values of m 
which are multiples of 4. In order to eliminate artifacts of the Trotter decomposition, 
we perform a binning of the histograms p(m) with four consecutive points in each 
bin. This implies that 

* = w < 3 - 12 > 

Altogether, we obtain 

—/ t \ p(m) M , . , m , 

= ~kr = Y p(m) - * 1 = 2M- (3 ' 13) 

with m constrained to be a multiple of 4. By construction, in the Euclidean time 
continuum limit the resulting probability distribution is normalized as 

d& x p($i) = 1. (3.14) 



Using the loop-cluster algorithm in its discrete-time variant jj, 14], we have 
simulated the spin ~ quantum XY model on a square lattice with L/a between 8 
and 64 at inverse temperatures (3 = L/c. We have worked at a sufficiently small 
lattice spacing in Euclidean time, such that the systematic discretization errors are 
negligible compared to the statistical errors. The probability distribution p($i) of 



11 




Figure 2: Probability distribution of the first component of the magnetization 

$1 on a 16 2 lattice obtained with the improved estimator. The error bars of the 
distribution are of the order of the line width in this figure. 
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the first component $1 of the magnetization has been obtained using the improved 
estimator described above. A typical distribution is shown in figure 2. As we will see 
below, the information about the vicinity of the minimum of the constraint effective 
potential is contained in the region of $1 where p($i) has its maxima. 

Due to the U(l) symmetry of the Hamiltonian, the probability distribution of 
the magnetization depends only on the magnitude of the order parameter $. 

Consequently, the probability distribution p($) can be expressed in terms of the 
probability distribution of the first component p(<&i) by using the relation 

/»27T /*0O 

p(<$> 1 )= dip d$ $ p($) 5($x - $ cosy?), (3.15) 
Jo Jo 

which can be cast into the form 



CO 



= 2 rf$ = p($). (3.16) 



($i/$) 



2 



The above relation is known as the Abel transform of the quantity p($) and can be 
inverted, provided that both p($) and its derivative p'($) tend to zero faster than 
1/$ as $ — > oo. For the probability distribution p($) one then obtains 

i r^}y / 1 (3.i7) 

Hence, given the Monte Carlo data for p($i), the probability distribution of the 
magnitude of the magnetization can be extracted. With the above equations one 
readily checks that the probability distribution p($) is properly normalized 

POO 

2tt / d$ $ = 1, (3.18) 
Jo 

provided that p($i) is normalized (see eq.f l3.14p ). We have determined the probabil- 
ity distributions p($) from p($i) using eq. (13.17jl . In figure 3 some representative re- 
sults for 27r$j9($) are shown. The mean value of $ decreases as the volume increases, 
and the width of the distribution p($) becomes narrower. It should be pointed out 
that the distribution is not symmetric around its maximum. The^ distribution turns 
into a 5-function in the infinite volume limit, centered at $ = Ai = 0.43561(1). 

Furthermore, we compute the first and second centered moments ($) and (($ — 
($)) 2 ) of the distribution p($). Once we have computed this can be done in 
a straightforward manner using the improved estimator which requires a computa- 
tional effort proportional to M 1,3 . 
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200 




Figure 3: Probability distributions 27r$p($) of the magnitude of the magnetization 
$ for L = 8a, 12a, 16a, 24a, and 64a. The lines are not fits but represent the Monte 
Carlo data themselves. The error bars are placed in equidistant positions o/$. The 
vertical line at $ = Ai = 0.43561(1) represents the ^-function distribution of the 
infinite system. 
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4 Comparison of Monte Carlo Simulations and 
Effective Theory Predictions 

In table 1 the first and second centered moments (<&) and (($— ($)) 2 ) of the distribu- 
tion are compared with the effective field theory predictions of eg. (12.251) . The 
errors of the theoretical predictions are due to the uncertainties in the low-energy 
parameters of eq. (12.311) and due to neglecting higher-order corrections. For the first 
moment the agreement is very good for L/a > 16. The absolute value of the second 
moment is very small and its statistical error is relatively large. Still, there are sys- 
tematic discrepancies between the Monte Carlo data of the first moment for small 
L, the second moment, and the 0(1/ L 2 ) effective theory predictions of eq. (j2.25j) . 
This discrepancy is well accounted for by additional 0(1/ L 3 ) corrections. Such cor- 
rections involve next-to-leading low-energy parameters which multiply higher-order 
terms in the effective action. At order 0(1/ L 3 ) we would also have to evaluate 3- 
loop graphs. This calculation, however, has not yet been worked out in the effective 
theory. Parameterizing the 3-loop terms with unknown coefficients a.\ and CK2, i.e. 



one obtains good fits to the Monte Carlo data for «i = —0.0013(2) and a 2 = 
—0.00061(5). This shows that the Monte Carlo data are described well by the theo- 
retical predictions. One may conclude that precise calculations of the two moments, 
combined with 3-loop predictions of the effective theory, would allow the determi- 
nation of some combination of sub-leading low-energy parameters. 

Starting from the probability distribution p($) one obtains the constraint effec- 
tive potential «($) by using the relation £>($) = A/"exp(— L 3 w($)). The constraint 
effective potentials corresponding to the curves in figure 3 are displayed in figure 
4. With increasing volume the constraint effective potential approaches the effec- 
tive potential, which is^ known to be a convex function. Using the rescaled variable 
ip = (pL / 'c)(<&— .M) / \M, one can also consider the extensive quantity U(ip) = L 3 u($) 
which is shown in figure 5. Expanding U(iJj) = U (tp) + (c/pL) Ui(ip) + 0(1/L 2 ), 
we have computed the universal part Uo(ip) by using Monte Carlo data for L/a 
between 8 and 64. Some values of the function Uo(ijj) extracted from the numer- 
ical data are compared with the analytic result of eq. (l2.17p in figure 6. It should 
be pointed out that the observed agreement does not rely on any adjustable pa- 
rameters. Even the normalization constant M of eq. (l2.13[) . which fixes an additive 
constant in the constraint effective potential, is predicted by the effective theory. 
As quantified in table 2, in the interval ip G [0,0.2], i.e. around the minimum of 
the constraint effective potential, the theoretical values of Uo(ip) and the numerical 
data agree remarkably well. To determine the low-energy constant /cq, we have used 




(4.1) 
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(($ - {<&)) ) thcory 


8 


0.46205(3) 


0.46224(1) 


1.90(7)xl0 4 


c\ o O 1 / 1 \ -i r\ — -4- 

2.7831(l)xl0 4 


12 


0.45305(10) 


0.45319(1) 


1.2(3) xlO- 4 


1.23694(6) xlO" 4 


16 


0.44875(3) 


0.44873(1) 


4.6(8) xlO" 5 


6.958(4) xlO" 5 


20 


0.44607(9) 


0.44608(1) 


4(2)xl0~ 5 


4.453(2) xl0~ 5 


24 


0.44432(10) 


0.44432(1) 


7(27) xl0~ 6 * 


3.092(2) xl0~ 5 


28 


0.44307(9) 


0.44306(1) 


4(22) xl0~ 6 * 


2.272(l)xl0~ 5 


32 


0.44198(8) 


0.44212(1) 


4(2) xl0~ 5 * 


1.7394(9) xl0~ 5 


40 


0.44086(9) 


0.44081(1) 


9(25) xl0~ 6 * 


1.1132(6)xl0" 5 


48 


0.43999(8) 


0.43994(1) 


6(200) xlO" 7 * 


7.731(4) xl0~ 6 


64 


0.43880(9) 


0.43885(1) 


4(250) xlO" 7 * 


4.3495(2) xl0~ 6 



Table 1: Comparison of Monte Carlo data (MC) for the first and second centered 
moments ($) and (($ — ($)) 2 ) o/p($) with predictions of the effective theory given 
by eg. 112. 25\) . The numerical errors of the analytical expressions are due to small 
uncertainties in the values of the low-energy parameters guoted in eg. Ii2.31\) . The 
entries with an asterisk (*) are statistically consistent with zero. The discrepancies 
between the Monte Carlo data and the effective field theory results are due to 3-loop 
corrections that were neglected in the theoretical predictions. 





#o(V0mc 


Uq{ip) theory 





- 2.985(6) 


- 2.980 


0.05 


- 3.781(4) 


- 3.781 


0.1 


- 4.229(2) 


- 4.232 


0.15 


- 4.222(2) 


- 4.224 


0.2 


- 3.632(5) 


- 3.640 



Table 2: Comparison of Monte Carlo data (MC) for the universal function Uo(ip) 
with the effective theory prediction of eg. \2.17\) . 
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Figure 4: Constraint effective potentials as functions of the magnitude of the 
magnetization $ for L = 8a, 12a, 16a, 24a ; and 64a. The constraint effective poten- 
tial approaches a convex effective potential in the infinite volume limit. 
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Figure 5: The extensive quantity U(<&) as a function of the rescaled variable ip = 
(pL/c)($ — for L = 8a, 12a, 16a, 24a, and 64a, compared to the analytic 

infinite volume result Uq(i[}). 
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0.05 0.1 0.15 0.2 



Figure 6: The analytic result for the universal function £/q(/0) compared with nu- 
merical values obtained from the Monte Carlo data for U(§). The comparison does 
not involve any adjustable parameters. 
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L/a 


B/J 


\{®){B)\ 


64 


0.00306135 


0.44053(3) 


64 


0.0042 


0.44142(4) 


72 


0.00328 


0.44076(4) 


72 


0.00386 


0.44121(7) 


80 


0.00266 


0.44028(4) 


80 


0.00313 


0.44069(5) 


88 


0.0022 


0.43985(3) 


88 


0.00258 


0.44021(6) 


96 


0.00184 


0.43950(3) 


96 


0.00217 


0.43984(5) 



Table 3: Monte Carlo data for \(Q)(B)\ which are used in the determination of ko- 

the extracted data for Ui(ip) as well as their theoretical prediction of eq. fl2.19p . A 
fit of the data to eq. fl2.19p leads to k = —0.0027(2). The result is illustrated in 
figure 7. Furthermore, ko can also be determined from fitting the -B-dependent field 
expectation values \(<&)(B)\ to their theoretical prediction of eq. fl2.24p . Since such 
an analysis provides a good check for the quantitative correctness of the numerical 
value for ko determined from Ui(ip), we calculate |($)(-B)| for several values of the 
magnetic field B and the box size L by using the technique of reweighting. Table 
3 contains the results of |($)(5)| obtained from reweighting. Using the data in 
table 3 as well as the corresponding theoretical prediction of eq.f l2.24l) . we arrive at 
k = —0.0026(3) which is in excellent agreement with k = —0.0027(2) calculated 
from Ui(ip). The statistical consistency between ko = —0.0026(3) determined from 
table 3 and the ko- value obtained from Ui(ip) also demonstrates the reliability of the 
reweighting technique employed for the determination of the data in table 3. Indeed, 
we have observed consistency between the data in table 3 with the largest 5-field 
for each L and the corresponding |($)(5)| determined by switching on explicitly a 
uniform magnetic field in the Monte Carlo simulations. Since for each L, the largest 
B imposes the greatest challenge for the reweighting method, we conclude that all 
the data in table 3 obtained by reweighting are indeed quantitatively correct. 



5 Conclusions 



We have computed the probability distribution of the magnetization in the (2 + l)-d 



XY model by using an improved estimator as first proposed in [13(] . Substantial 



improvements in the implementation of the method have allowed us to reach spatial 
volumes as large as 64 2 (or even 96 2 to calculate \(§)(B)\). Using the improved 
estimator in a loop-cluster algorithm simulation, we have determined the first and 
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Figure 7: Result of fitting the Monte Carlo data for Ui(ip) to the theoretical predic- 
tion. 
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second moments (<£>) and (($ — ($)) 2 ) of the distribution p(<£>) of the magnitude $ 

— # 

(of the mean magnetization vector $), as well as the constraint effective potential 
■u($) (obtained fromp(<3>) = A/"exp[— L 3 -u(<3>)]) for different space-time volumes. The 
Monte Carlo data are in excellent quantitative agreement with analytic predictions 
which Gockeler and Leutwyler derived from a systematic low-energy effective field 
theory. This shows that the magnon effective field theory indeed provides a quanti- 
tatively correct systematic derivative expansion of the low-energy physics. Thanks 
to the very efficient loop-cluster algorithm, in the context of the (2 + l)-d quantum 
XY model we were able to test theoretical predictions of the effective field theory 
up to two-loop order. 
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